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1. Introduction 
1.1. Characteristics of microalgae 


With the continuous reduction of the fossil fuels, and the 
acceleration of global greenhouse effect mainly due to carbon 
dioxide (CO2) emissions, energy supply may be in trouble in the 
near future. It is urgent and significant to find reliable clean energy 
resources alternative to fossil fuels [1]. Microalgae may become 
one of the most promising new resources to supply energy and 
mitigate CO2 in the future [2,3]. Technologies for producing 
microalgae and using microalgae for biodiesel have been known 
for more than 50 years. Up to now, many researchers have done 
plenty of work on planktons [4-7]. Compared to the first genera- 
tion biofuels, microalgae have several advantages in sustainability, 
economics and environment. Microalgae not only have higher 
productivity, but also can be fed in saline/brackish water/coastal 
seawater on non-arable and deserted land. However, there are 
some shortcomings of microalgae for use as biofuels. For example, 
the dry solids content of the microalgae and water mixture can be 
as low as 0.05% dry solids content. The energy balance associated 
with segregating the water from the microalgae dry solids to allow 
a bio-esterification process to be undertaken is very high [8]. 
Literature [9] also indicated that some unstable biodiesel with 
many polyunsaturates will be derived in the process of producing 
biodiesels from algae oil. 


1.2. Microalgae models 


Some models have been established by various researchers 
to predict different aspects of microalgae growth process, such 
as mass growth, light distribution, pigment dynamics, nitrogen 
uptake, depth dependence, temperature dependence and respira- 
tion rate, etc. [10-13]. The biological model of microalgae comes 
from the population model proposed by Malthus [14]. The concept 
of the population model such as growth rate is used in later 
microalgae growth models. Droop [15,16] proposed a dynamic 
model which takes the dilution rate and influent inorganic nitro- 
gen concentration into account to describe the growth of micro- 
algae. The model is classical, but practical, which describes the 
kinetic model of the algae growth and the nitrogen uptake. 
Recently, several new models were proposed based on Droop's 
model. Geider et al. [17] proposed a new model which included 
growth process, nitrogen uptake, chlorophyll synthesis, tempera- 
ture and respiration aspects, but depth dependence was ignored. 
Quinn et al. [18] proposed a model for industrial scale systems 
based on Geider et al.'s model [17]. A dynamic model was 
proposed by Packer et al. [19] to predict the growth and neutral 
lipid synthesis of green algae by using a new method. This method 
considered the influences of two factors (i.e., the photosynthesis, 
and the nitrogen uptake) on the growth rate, which was defined as 
a piecewise function of the two factors. Bernard [20] proposed a 
different kinetic model which was verified by experimental data 
[21]. Different from Packer et al.'s model [19], the light-limited and 
nitrogen-limited factors will simultaneously have effects on the 
growth rate in the model. Mairet et al. [22] further developed the 
previous models by dividing the biomass into three compartments 
including the function part, carbohydrates, and neutral lipids. The 
transformations and balances between the three parts were dis- 
cussed in their model. 


In the above models, various expressions and coefficients have 
been used. Investigation on the differences between them is needed 
so that it is easier for later researchers to select more proper 
expressions and coefficients to model, predict and control the 
microalgae growth. By now, little work on the comparison of 
different microalgae sub-models (that is, the expressions or coeffi- 
cients, for example, the sub-model of nitrogen uptake) has been 
reported. In this paper, based on the valid Bernard's dynamic model 
[20] acting as the basic model (not the standard model) made up of 
several sub-models, various expressions and coefficients in several 
main models [17-19] will be compared by comparing the basic 
model and a new model. The new model is produced when one 
sub-model of the basic model is replaced by the corresponding 
sub-model used in another model. In the second and the third 
parts, the influence of other aspects such as the depth of culture and 
the ambient temperature is also discussed. 


2. Model description 


Many models have been proposed up to now, which used 
various expressions to calculate the parameters e.g. light distribu- 
tion, pigment dynamics, nitrogen uptake, growth rate, respiration 
rate, and temperature dependence. In this part, we will describe 
Bernard's model, and the various expressions of parameters for 
various models [17-20] are also discussed. 


2.1. Bernard's model 


The model [20] describes four variables in the ordinary diffe- 
rential equations: s(t), which denotes the concentration of dis- 
solved inorganic nitrogen, the nitrate or ammonium (g N m7), 
q(t), which is internal nitrogen cell quota (g N (g C)~!), x(t) which 
is algae biomass concentration (g C m~3), and [*(t) which is not 
the real radiation but a conceptual variable denoted radiation. The 
unit of [*(t) is mol m~? s~!, which means the number of photons 
absorbed per unit area per unit time. The expression of calculating 
variable I* will be discussed later. The four ordinary differential 
equations were expressed as: 


A ian 5 q 

S = DSin rkl E) Ds (1) 
š _ S q = m 

i= a) (Io, 1", x, 4)(q — Qo) (2) 
„=g pk Qo 

X= Allo, I*, x, q) Ea x—Dx—Rx (3) 
P =F, I*,x,q) (1 - 2) a-r") (4) 


where D denotes the dilution rate, p is the maximum nitrogen 
uptake rate, Sin is influent inorganic nitrogen concentration, R is 
the inspiration rate, J is the average radiation along culture 
volume, Qj is the minimum nitrogen quota, Q4 is the maximum 
nitrogen quota, and z denotes the average growth rate, which is 


calculated by 


2K Ip(1 —e7*)V/A 


Ho, é) = À il arctan 
g AVA 2Re- +o +e-4)Kiy + 2225.95) 


(5) 
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where lo represents the light intensity on the bioreactor surface, 
and ñ is the maximum growth rate. The average radiation along 
culture volume I is given by 


=Í 
=70 -e7 (6) 


where the optical depth, 2, is the product of the depth of the 
culture L multiplied by light attenuation rate é 


ASL (7) 


The light attenuation rate é can be calculated by 
&€=aChl+bx+c (8) 


where, a, b, and c are the constants, and the chlorophyll concen- 
tration Chl can be calculated by 


Chl = y(I")xq (9) 


The proportion of chlorophyll concentration to nitrogen con- 


centration y(I*) can be given by 
MO) = tare (10) 
mane Ke 


where ymax is the maximum value of y(I*). 
The parameters A and I,,,which denote the radiation providing 
maximum rate of photosynthesis in Eq. (5) are calculated by 


A =4Ê 


opt 


lope = Ks Kiy (12) 


where Ks is 


(0S) —K?, (11) 


Ky = =! (13) 


The initial value of Chl/x is denoted by a parameter 6§ 


Chlo 


iC 
05 = 
0 Xo 


(14) 


where Chip is the initial value of Chl, and xg is the initial value of x. 


2.2. Light distribution modelling 


Light intensity is an important factor of the algae growth model 
[23]. In nearly all the models, the Lambert-Beer Law was used to 
describe the light distribution in the reactor. The law assumes that 
the light below the culture surface is attenuated at an exponen- 
tially decaying rate because the biomass and/or the chlorophyll 
absorb a part of the radiation from the culture surface 


KD) = Ipe~ (15) 


Eq. (6) can be obtained from Eq. (15). The exponential factor of 
light attenuation in the culture, é, may be related to the biomass 
content and/or the chlorophyll content. There are four different 
expressions of calculating £: (1) a # 0, b + 0, c = 0, that is, the factor 
é depends on both the biomass concentration and the chlorophyll 
concentration, as reported by [20]; (2) a=b=0, c#0, that is, 
parameter é=c when the reactor depth L is not too deep, as 
reported by [17]; (3) b #0, a= c= 0, that is, the light attenuation 
rate linearly correlated to biomass concentration, as reported by 
[18,24]; and (4) a # 0, b = c = 0, that is, light attenuation coefficient 
é is linearly related to chlorophyll content, as reported by [11,19]. 


2.3. Pigment dynamics modelling 


Besides the known a, b, and c constants, the chlorophyll 
concentration in the reactor must be known for calculating the 
light attenuation coefficient. Bernard [20] assumed that the 
chlorophyll concentration is proportional to the nitrogen con- 
centration, xxq, as demonstrated by Eq. (9). The proportional 
coefficient is calculated with Eq. (10). Quinn et al. [18] considered 
a=O in Eq. (8). It is not necessary to calculate the chlorophyll 
content for calculating the light attenuation coefficient. Several 
other methods were also used to calculate the chlorophyll 
content. 


(1) From Geider et al.'s model 
In this model [17], the chlorophyll content is calculated with 
the following differential equation 


1 dChl _ hcnp 
Chl dt o 
where p is the uptake rate of nitrogen, @°is the ratio of 
chlorophyll a to phytoplankton carbon, rcp, is the chlorophyll 
a degradation rate constant, and hc, is the chlorophyll a 
synthesis regulation term, which is calculated by 


—Tchi (16) 


x 


hen = N = 17 
Chl max Chl gC To (@ ) 


where 8" is the ratio of chlorophyll a to N, 6%,,, is the maximum 
value of 6", and a“! is the chlorophyll a—specific initial slope of 
the photosynthesis-light curve. 

From Packer et al.’s model 

The model [19] used the chlorophyll quota 6° to calculate the 
chlorophyll content 


~ 
N 
i 


of = Core Or (18) 


where Co is the carbon subsistence quota, and p is the photo- 
synthesis rate. Eq. (18) didn't take the respiration into 
account. 


2.4. Nitrogen uptake modelling 


Nitrogen is an important material in microalgae culture. The 
nitrogen source in bioreactors can be inorganic source such as 
ammonium or nitrate and organic source. Microalgae can synthe- 
size proteins and neutral lipids by absorbed nitrogen. It is 
necessary to model the nitrogen uptake. In Bernard's model [20] 
and Packer et al.'s model [19], the nitrogen uptake behavior was 
described with Eq. (2). Then a brief introduction to the origin of 
Eq. (2) will be given. 

The nitrogen cell quota of the culture, q, is equal to the ratio of 
the nitrogen concentration to the carbon concentration. This is 
analogous to the C:N ratio which is of crucial importance to the 
digestibility of substrates. q can be expressed as: 

S . Sx-xs Š X 
1=4% 74 x2 x x 


p(s, q) — HU. q, £q (19) 


where s=xxq, and yw is the growth rate. q changes due to the 
nitrogen uptake and the microalgae growth. The right side of Eq. 
(19) can be divided into two parts: the nitrogen uptake part p(s, q) 
and the growth dilution part w(o, q, 8&q 


8.9) =Ps¢K.(1-9.) (20) 
lo. 4. |= Hlo.2(1 -22) (21) 
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where p is the maximum nitrogen absorption rate. Substituting 
Eqs. (20) and (21) into Eq. (19), Eq. (2) is derived. Then several 
methods to model the nitrogen uptake are shown: 


(1) From Geider et al.’s model 
In this model [17], Nitrogen uptake is described with a diffe- 
rential equation. 


>P 
P 22 
q q Ty (22) 
ae Qi-4q 

PEPEK: G rm (23) 


where ry is the remineralization rate constant, and f(T) 
denotes the temperature function. 

From Quinn et al.’s model 

In this model [18], the temperature function f(T) has been 
considered. The equation of nitrogen uptake in Quinn et al.’s 


= 
N 
paia 


model is 
gp T 24) 
S+K, q 


where Ry is the respiration constant for nitrogen. 


2.5. Growth rate modelling 


It is universally acknowledged that the growth rate is correlat- 
ing to the radiation intensity. In Bernard's model, the nitrogen- 
limited and the light-limited are two multiplication factors, as 
shown by Eqs. (5) and (21). Eqs. (5) and (11)-(14) are used to 
calculate the growth rate of microalgae. There were some methods 
to calculate the microalgae growth rate in various models. Some 
models considered the light radiation effects on photosynthesis 
rate, others (e.g. Packer et al.’s model [19] and Bernard's model 
[20]) related the growth rate to the radiation intensity directly. 
In Packer et al.’s model [19], the growth rate was determined by 
both photosynthesis rate and the nitrogen uptake rate. 


(1) From Geider et al.’s model 
Geider et al. [17] assumed that the growth rate is the result of 
subtracting the energy consumption from the photosynthesis 


1d 
== p—Re-p (25) 
— hl gc 
p= Pr exo( 2 An) (26) 
Pm 
pE i) 27 
Pm Phy (a ( ) 


where Rc is the maintenance respiration rate, ¢ is the cost of 
biosynthesis, p,, is maximum value of p at temperature T, and 
PS is maximum photosynthesis rate. This model assumed that 
50% of the biomass is carbon. 

(2) From Quinn et al.’s model 
Quinn et al.’s model [18] is the same as the Geider et al.’s one 
[17] except the use of different method to calculate pm 


Pm = Pref (1 -2%) (28) 


The nitrogen uptake factors expressions in Eqs. (27) and (28) 
are different. After comparing the variations of the values of 
(q—Qo0)/(Q1 — Qo) term in Eq. (27) and the values of (1 — Qo/q) 
term in Eq. (28) with nitrogen quota q varying from 0.05 g N 
(gC)~! to 0.25gN(gC)~|, it is found out that the nitrogen 


uptake factor in Eq. (27) linearly increases from 0 to 1, and that 
in Eq. (28) logarithmically increases from 0 to 0.8. 

From Packer et al.’s model 

The growth rate is interestingly shown as a piecewise function 
of time. The microalgae growth rate is determined by both 
nitrogen and light factors. The growth rate is therefore 
considered to be the minimum of the growth rates calculated 
by using the nitrogen-limited and light-limited relationships 


~ 
Ww 
= 


[19]. 
X = ux — Dx (29) 
u= minfa (1-2).2} (30) 
p=OP, Í exp (+ “| (31) 
(qx) pS 
m= MEL (32) 
(qx)* +(Qox) 


where the value of pC, is selected to be 7 day~ 1 in Eq. (32), A is the 
optical cross section of chlorophyll a, Cp is the carbon subsistence 
quota, ji is 3.26 day~' in Eq. (30), and @ is the Quantum efficiency. 


2.6. Respiration rate and degradation rate modelling 


There are five sub-models to model the respiration or the 
degradation of microalgae. In Bernard's model [20], respiration is 
taken into account only in the growth modelling. Quinn et al. [18] 
proposed a complex method to distinguishing two forms of respira- 
tion: respiration of growth, Rc, and the reapiration of nitrogen uptake, 
Ry. The two parameters Ry and Rc are used in Eqs. (24) and (25) to 
represent the two forms of respiration. 

In Geider et al.’s model [17], three respiration or degradation 
rates, chlorophyll a degradation rate constant, rcp), nitrogen 
remineralization rate constant, ry, and Rc were used. The three 
rates are equivalent, Re=ry =rcn =0.025 day~!, which are used in 
Eqs. (16), (22) and (25), respectively. 

The respiration or degradation rate is considered as a constant 
in all the above sub-models. Mairet et al. [22] didn't regard it as a 
constant respiration rate, but defined the respiration rate propor- 
tion by using a linear relationship to the nitrogen uptake rate p 


R=19+11p(S,Q) (33) 


where rg and r; are the coefficients. 


2.7. Depth dependence modelling 


In massive production, depth is a vital parameter in bioreactors 
such as in raceway or pond. If the bioreactors are deeper, more 
light radiation can be available for use, but more materials will be 
wasted. The influence of the variable depth will be discussed in 
next part. 


2.8. Temperature dependence modelling 


The function related to temperature T consists of two equations 
which were proposed by Quinn et al. [25]. 


26, 

D= 34 

fT d+) (34) 
Ea _ Ea 

Dy = e Fo! (35) 


where, Ro is the universal gas constant, Top is the optimum 
microalgae growth temperature, and E, is the activation energy 
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carboxylation Rubisco. Now the remaining problem is how to add 
the temperature function to the ordinary differential functions. 
The temperature function can be added to Bernard's model [20] by 
assuming the nitrogen uptake rate and the average growth rate to 
be related to the temperature factor. The equations of modified 
model of Bernard's are 


; Ss q 

$= DSin Pet A) Dx Ds (36) 
a=r g (1-0 JAD Allo. x. Da -fT 67) 
X = fill, I*, x, o(1 - So) fx- Dx- RFD (38) 


When the environmental temperature T is equal to the opti- 
mum temperature Topt, the temperature function f(T)=1, thus the 
model will degenerate to Bernard's model. 

In Quinn et al.'s model [18], the temperature factor is taken into 
consideration. Substituting temperature function into Eqs. (20), 
(27) and (28), the below equations are derived 


->s fy 8. 
8.9) =P l -EVO G9 
c q— Qo ) 
Pm Pref G =Qg fT) (40) 
Pm = Pref (1 — ©) Far) (41) 


Based on Bernard's model [20], a modified Quinn et al.’s model 
[18] is established whose growth part and temperature function are 
substituted by Eqs. (39)-(41) in Quinn et al.'s model. The model will 
be discussed in the next part. 

Geider et al. [17] also took temperature effect into consideration 
in nitrogen uptake, as shown in Eq. (23). Packer et al. [19] didn’t 
take the temperature factor into consideration in their model. 


S 
0.15 


— Simulation 
o Experiment} | 


0.1: eny 


NH4(mol/m?) 
fo) 
O 
o 


10 20 30 


Chl 


Chl(g Chi/m’) 


30 
Time(d) 
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3. Results and discussion 


Fig. 1 compares the simulations using this model and the 
experimental data about Isochrysis galbana performed by Flynn 
et al. [21]. The parameters used in the simulations with the model 
are presented in Table 1. It is assumed that no additional nutrition 
is added into the culture once the growth of microalgae began. 
As shown in Fig. 1, the model simulations have a good agreement 
with the experiment data [21], which shows the validation of the 
simulations. This figure indicated that the extracellular nitrogen 
was exhausted on the 12th day. At the same time, the synthesis of 
chlorophyll ceased. The growth process of biomass content with 
time could be divided into four phases [21]. The first five days 
were the lag phase, during which the cells increased at a slow rate. 
Following the lag phase was the exponential phase (from the 5th 
day to the 12th day), during which the growth rate of cells was 
0.4 day~'. The cells didn’t stop growing until the 19th day, and 
these days (from the 12th day to the 19th day) were the post- 
exponential phase. After that, the culture entered the stationary 
phase, during which the content of biomass, chlorophyll and 
intracellular nitrogen were stable. 


3.1. Light distribution modelling 


In order to illuminate the influence of various light distribution 
relationships on the simulation results, simulations based on 
Bernard’s model by using different expressions of the light 
attenuation coefficient € are performed and compared. 

The simulations show that, for L=0.05 m, compared with light 
distribution sub-models proposed by Quinn et al.’s [18] and Packer 
et al.’s model [19], the chlorophyll values for sub-models from 
Bernard’s model [20] and Geider et al.’s model [17] are a bit lower 
in algae biomass concentration while the inorganic nitrogen 
concentration, chlorophyll concentration and intracellular nitro- 
gen concentration are almost the same. Things are different when 
the reactor becomes deeper. When L=0.5 m, the differences of the 
four curves corresponding to the four expressions of the light 
attenuation coefficient € become obvious though their variation 


30 


0 4 
0 10 20 30 


Time(d) 


x*q 


0 10 20 30 
Time(d) 


Fig. 1. Comparison of simulation and experimental data from [21]. 


530 


trends with time are the same. Compared with light distribution 
sub-models proposed by Quinn et al.’s [18] and Packer et al.’s 
model [19], the growth rate and nitrogen uptake rate results for 
sub-models from Bernard’s model [20] and Geider et al.’s model 
[17] are lower while the chlorophyll content is even larger. 


3.2. Pigment dynamics modelling 


Simulations based on Bernard’s model [20] by using different 
expressions of the chlorophyll content are performed and com- 
pared. Fig. 2 shows that the variations of four main parameter 
values with time calculated based on Bernard’s model by using 
three different expressions of the chlorophyll concentration Chl 
when L=0.05 m. 


Table 1 
Parameters values used in simulation. 


Parameters Values Units 
ji 17 Day~! 
Qo 0.05 gN (gc)! 
Qi 0.25 gN(gC)"! 
Ki 14 mol m~?s~! 
Ki 295 mol m~? s~! 
7 0.073 gN (gC)! day™! 
K, 0.0012 gNm? 
R 0.0081 Day~? 
Ymax 0.68 g Chl (g N)`! 
kp 63 mol m~? s~! 
ach! 0.46 m? (g Chl)~! 
Ip 100 molm~?s~! 
4 2.0 gC(gN)"! 
ON ox 0.3 gChl(gN)7! 
Poy 3.0 Day"! 
Co 0.61 gC(gdw)7! 
A 4.82 m? g~! Chl 
D 9.84 x 107? g C (mol photons)~! 
Ea 63 kJ mol`! 
Ro 8.3145 JK! mol"! 
Toe 18 fa 
S 
0:15 


— Eqs. (9) and (10) 
——Egs. (16) and (17) 


0 10 20 30 
Time(d) 


Chl 


Chl(g Chi/m’) 


30 
Time(d) 
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As shown in Fig. 2, the dissolved nitrogen concentration, and 
the carbon and nitrogen contents for different pigment sub- 
models are almost the same expect chlorophyll content. Fig. 3 
shows the four main parameter values when L changes to 0.5 m. 
The differences in inorganic nitrogen concentration, intracellular 
carbon concentration and intracellular nitrogen concentration 
results of different pigment sub-models appear and the differences 
between the chlorophyll content are still evident. 


3.3. Nitrogen uptake modelling 


Fig. 4 shows that the variations of four main parameter values 
which is same as Fig. 2 with time calculated based on Bernard’s 
model by using three different expressions of the nitrogen uptake 
when L=0.05 m. The chlorophyll content, intracellular carbon 
concentration and intracellular nitrogen concentration result for 
nitrogen uptake sub-models from Geider et al.’s model [17] is 
smaller than the sub-models from Bernard’s model [20] and Quinn 
et al.’s model [18] while the inorganic nitrogen concentration is 
almost the same. 


3.4. Growth rate modelling 


Fig. 5 shows that the variations of four main parameter values 
which is same as Fig. 2 with time calculated based on Bernard’s 
model by using four different expressions of the growth rate when 
L=0.05 m. Note that parameter Prep is selected as 35 day~! in 
Geider et al.’s model and 30 day~! in Quinn et al.’s model. 

As seen from Fig. 5, it is concluded that the four curves have the 
same trend. The growth rate for growth sub-models from Geider 
et al.'s model [17] and Quinn et al.’s model [18] are drastically 
larger than the sub-models from Bernard’s model [20] and Packer 
et al.’s model [19] after the fourth day. But the inorganic nitrogen 
concentration, chlorophyll content and intracellular nitrogen con- 
centration result for growth sub-models from four references 
approximate to each other. 
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Fig. 2. Three different expressions to calculate the chlorophyll content when L=0.05 m. 
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3.5. Respiration rate modelling 


Fig. 6 shows that the variations of four main parameter values 
calculated based on Bernard’s model [20] by using five different 
expressions of the respiration rate when L=0.05m. Based on 
Bernard’s model, the respiration parts in every equations are 
added in Eqs. (2) and (3). The five cases are analyzed as follows: 
(1) adding — Rx into Eq. (3), where R=0.0081 day~'; (2) assuming 
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R=0; (3) adding —Rcx into Eq. (3), where R-=4.32x10~4 day~!; 
(4) adding — Rx in Eqs. (2) and (3), where R=0.0081 day~!; and 
(5) calculating the factor R with Eq. (33). 

As shown in Fig. 6, the growth rate for respiration sub-model 
proposed by Mairet et al. [22] is a bit smaller than the sub-models 
from the other four. In general, the respiration factor is not an 
important aspect in the whole model, but an accuracy method is 
necessary if we want to get a better result. 
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Fig. 3. Three different expressions to calculate the chlorophyll content when L=0.5 m. 
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Fig. 4. Four main parameter values calculated based on Bernard's model by using three expressions of nitrogen uptake when L=0.05 m. 
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3.6. Depth dependence modelling 


In this paper, the depths of bioreactor, L, are 0.05 m and 0.5 m, 
which represent the shallow bioreactors and the deep bioreactors, 
respectively. The variations of the four parametric values with 
different depths are discussed and compared in Fig. 7. As shown in 
Fig. 7, the results of four main parametric values vary a little with 
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the depth when the depth L varies from 0.03 m to 0.15 m. This 
shows that it is reasonable to take L=0.05m as a reference 
example to represent the shallow depth in this paper. Fig. 8 shows 
four parametric values vary with time when the bioreactors 
become deeper. The results of four main values vary a lot with 
the depth when the depth L varies from 0.05 m to 1.00 m. When 
L=1.00 m, not only the biomass concentration become lower, but 
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Fig. 5. Four main parameter values calculated based on Bernard’s model by using four different expressions of growth rate when L=0.05 m. 
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Fig. 6. Four main parameter values calculated based on Bernard’s model by using five different expressions of the respiration rate when L=0.05 m. The five expressions are: 
(1): adding — Rx into Eq. (3), where R=0.0081 day—!; (2): assuming R = 0; (3): adding —Rcx into Eq. (3), where Re =4.3210~* day !; (4): adding —Rx in Eqs. (2) and (3), 
where R=0.0081 day~!; (5): calculating the factor R with Eq. (33). 
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also more materials such as nitrate are demanded. Compared 
these two figures, it is found out that the biomass concentration 
will decrease a lot when the PBR depth is larger than 0.15 m. 


3.7. Temperature dependence modelling 


Fig. 9 shows that the variations of four main parametric values 
calculated based on the modified Bernard’s model [20] by using 
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five different temperature rates when L=0.05m. As shown in 
Fig. 9, based on the assumption of the optimal temperature of 
18°C for Isochrysis galbana, the modified Bernard’s model is 
degenerated to Bernard’s model without considering the tempera- 
ture factor. Although there are same trends for all the five curves at 
various temperatures (T=8, 16, 18, 20, and 28 °C), microalgae grow 
slower when the environment temperature is deviated from the 
optical temperature of the algae. 
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Fig. 7. Four main parameter values calculated based on Bernard’s model by using four different expressions at five depths (L=0.03 m, 0.05 m, 0.08 m, 0.10 m, 0.15 m). 
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Fig. 8. Four main parameter values calculated based on Bernard’s model by using four different expressions at five depths (L=0.05 m, 0.10 m, 0.30 m, 0.50 m, and 1.00 m). 
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Fig. 9. Four main parameteric values calculated based on modified Bernard’s model at five different temperatures. 


4. Conclusions 


Microalgae acted as the second generation biofuels could be 
produced in quantity with culture technologies in artificial bior- 
eactors. Modelling microalgae growth is an effective method to 
predict the growth process, and then estimate the production and 
help to improve the productivity of microalgae under a given 
culturing condition at a given site. Recently, some models have 
been proposed with various expressions and coefficients of para- 
meters. In this paper, the differences due to use of the various 
expressions and coefficients of parameters in the models have 
been compared. The conclusions are drawn (based on the case 
L=0.05 m) as follows: 


(1) Compared with light distribution sub-models proposed by 
Quinn et al.’s [18] and Packer et al.’s model [19], the chlor- 
ophyll values for sub-models from Bernard’s model [20] and 
Geider et al.’s model [17] are a bit lower in algae biomass 
concentration. The inorganic nitrogen concentration, intracel- 
lular carbon concentration and intracellular nitrogen concen- 
tration results of different pigment sub-models are almost the 
same expect chlorophyll content. 

(2) The chlorophyll concentration, intracellular carbon concentra- 
tion and intracellular nitrogen concentration result for nitro- 
gen uptake sub-models from Geider et al.’s model [17] is 
smaller than the sub-models from Bernard’s model [20] and 
Quinn et al.’s model [18] while the inorganic nitrogen con- 
centration is almost the same. 

(3) The growth rate for growth sub-models from Geider et al.’s 
model [17] and Quinn et al.’s model [18] are drastically larger 
than the sub-models from Bernard’s model [20] and Packer et 
al.’s model [19] after the fourth day. The growth rate for 
respiration sub-model proposed by Mairet et al. [22] is a bit 
smaller than the sub-models from the other four. 

It is concluded that the biomass concentration will decrease a 

lot when the PBR depth is larger than 0.15 m. Microalgae grow 

slower when the environment temperature is deviated from 
the optimum temperature of the algae. 


S 


This paper will guide the researchers to select the most 
appropriate models or sub-models or build new models in order 
to model the microalgae growth process. 


Nomenclature 

A optical cross section of chlorophyll a (m? g~! ChI) 

Co carbon subsistence quota (g C (g dw)~') 

Chl chlorophyll concentration (g Chl (g dw)~!) 

D dilution rate (day~ ') 

Eq activation energy carboxylation Rubisco (J mol~') 

f(T) temperature function (Dimensionless) 

F nitrogen restriction factor (Dimensionless) 

hcni chlorophyll a synthesis regulation term (Dimensionless) 

Io light intensity at reactor surface (mol m~? s7!) 

Iopt radiation providing maximal rate of photosynthesis 
(mol m~? s~!) 

O conceptual radiation (mol m~? s7!) 

I average radiation along culture volume (mol m~? s~!) 

L depth of culture (m) 

p photosynthesis rate (day~') 

Dis maximum value of p at temperature T (day~!) 

Prep maximum photosynthesis rate (day~ ') 

q(t) internal nitrogen cell quota (g N (g C)~!) 

Qo minimum nitrogen quota (g N (g C)~!) 

Q1 maximum nitrogen quota (g N(gC)~') 

R respiration rate (day~!) 

Ro universal gas constant (J K7! mol~!) 

Rc maintenance respiration rate (day~') 

Tchi chlorophyll a degradation rate constant (day~!) 

Ty nitrogen remineralization rate constant (day~!) 

Ry respiration constant for nitrogen (day~ ') 


concentration of dissolved inorganic nitrogen, (nitrate or 
ammonium) (g N m~?) 

Sin influent inorganic nitrogen concentration (g N m~?) 
optimum growth temperature (°C) 

biomass concentration (g C m~?) 
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ach! chlorophyll a-specific initial slope of photosynthesis- 
light curve (m? (g Chl)~') 
yl") proportional coefficient of Chlorophyll concentration: 


nitrogen concentration (g ChI(g N)~') 
Ymax maximum value of y(I*) (g Chl (g N)~') 


¢ cost of biosynthesis (g C (g N)~!) 

0s initial value of 6° (g Chl (g C)~!) 

oF chlorophyll a: phytoplankton carbon ratio (g Chl (g C)~') 
oN maximum chlorophyll a: N (g Chl (g N)~') 

ON ox maximum value of 6% (g Chl (g N)~') 

A optical depth (Dimensionless) 


u growth rate (day~') 
m average growth rate (day~') 

it maximal growth rate (day~ !) 

€E light attenuation rate (m~!) 

p uptake rate of nitrogen (g N(gC)~! day~!) 

p maximum nitrogen absorption rate (g N (g C)~! day~ !) 
p Quantum efficiency (g C (mol photons)~')) 
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